home *** CD-ROM | disk | FTP | other *** search
Text File | 1991-02-19 | 469 b | 27 lines | [TEXT/????] |
- Sx = 0
- Sy = 0
- Sxy = 0
- Sx2 = 0
- Sy2 = 0
- n = 5
- I = 1
- DO WHILE(I<6)
- Sx = CR(1,I) + Sx
- Sy = CR(2,I) + Sy
- Sxy = CR(1,I) * CR(2,I) + Sxy
- Sx2 = CR(1,I)^2 + Sx2
- Sy2 = CR(2,I)^2 + Sy2
- I = I + 1
- REPEAT
- TEMP = n*Sxy - Sx*Sy
- B = TEMP / (n*Sx2 - Sx^2)
- A = (Sy - B*Sx)/n
- r = TEMP / SQRT((n*Sx2 - Sx^2)*(n*Sy2 - Sy^2))
- AREA1 = %DFIT(1,2,5)
- AREA2 = %LFIT(A,B,CR(1,1),CR(1,5))
- IF (AREA2 > AREA1)
- PError = 100.0*(AREA2 - AREA1)/AREA1
- ELSE
- PError = 100.0*(AREA1 - AREA2)/AREA2
- END
-